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ABSTRACT 

The rate of breakup of helium nuclei by particle-induced 
reactions is computed. It is shown that the rate is 
determined by the endothermic reaction p + He He + d 
beco mi ng effective at kT > few MeV. It is suggested 
that matter having been processed to these temperatures 
will be depleted in helium and in the elements C, N, 0, 
and Ne. 


1. Introduction. In recent years considerable interest has developed 
in relativistic plasmas in which the electron gas is highly relativis 
tic (0 = kT/mc 2 * 10). At these temperatures the ions would be non 
relativistic but capable of undergoing endothermic reactions. In a 
recent paper [1] devoted to nuclear-force effects in p-p scattering, it 
was suggested that helium nuclei will break down in the following series 
of endothermic reactions: 

p + 4 He + 18.353 MeV -> 3 He + d (1) 

p + 3 He + 5.493 MeV +d+p+p ^ 

p + d + 2.225 MeV -»• p + p + n 


Because of its significantly larger threshold, the reaction (1) would be 
the slowest and therefore determine the rate of the overall helium de 
struction process. Once the helium is broken down to hydrogen, any sub- 
sequent buildup would be very slow because of the small cross section 
M ^ J 1 e + v rt . The plasma would 


d + 


e‘ 


for the fundamental p-p reaction p + P c , , . 

then remain essentially pure hydrogen. Other elements would also be 
broken down (see [1]), such as the abundant medium-weight C, N, 0, Ne, 
although large Coulomb barriers are involved, since the threshold ener 
gies for the breakup reactions are much larger than in reaction (I) . 


The most elementary estimate of the characteristic breakdown 
temperature is obtained [1] by setting kT 'v 0 >1E b > E b 

binding energy difference (18,353 MeV) in reaction (1). The factor 0.1 
arises since the equation for the breakdown temperature will be of the 
form, characteristic of ionization and dissociation phenomena, 

E. /kT //•> 

e N = large number . v / 


The logarithm of any large number is ^ 10, and this yields for the 
breakdown temperature: 

0 b (He) ^ 4 ( 

2. Calculation of Breakup Reaction. The calculation of the character 
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istic breakdown temperature is really a non-equilibrium problem requir- 
ing a knowledge of the density (n ) of the plasma and the available 
time (t) . If r = dN/dVdt is the number of reaction (1) per unit 
volume per unit time, the breakdown temperature will be computed from 


The rate r will be proportional to n 2 and an integral over the reac- 
tion cross section: ” 



where v r is the relative velocity of the reacting protons. 


The problem is in the correct use of the reaction cross section 
o and in the extrapolation of its experimental values to low energies. 
For in this case of a reaction with a threshold energy and reactants 
with a Maxwellian velocity distribution, the main contribution to 
< crv r > will come from c.m. energies just above threshold. However, 
then the Coulomb effects in the cross section will be very important, 
especially for the outgoing particles in the reaction. To include these 
effects the experimental cross section [2] has been fit to a form 


where CT nuc is a nuclear part" and C 2 and C' 2 are (incoming and 
outgoing) Coulomb factors of the form [3] 

C 2 = 2TTn/(e 2lTT1 -l) ( 

for the incoming and outgoing particles. Here n = z,z ? (Ry /E) 1 / 2 , in 
terms of the charges z^e and Z£e of the particles involved and 
R y y = (y/m) Ry is the Rydberg energy corresponding to the reduced mass 
of the pair of nuclei. In terms of the c.m. energy (before reacting) 

and the threshold energy E fc for reaction (1) , a is found to have 
the form nuc 


°nuc = °nuc-t [1 " ^"V^t 1 (10) 

wit h o nuc _ t = 93mb and a g = 0.148. 

Due primarily to the Coulomb factor C' 2 for the outgoing 
particles, the integrand in < ov r > is found to have a "Gamow peak" 
where most of the contribution is centered. An asymptotic evaluation 
of the integral yields the following set of breakdown temperatures for 
given values of xn p : 

0 b (He) 2345678 

xn p (yr cm" 3 ) 1.4xl0 15 3.8xl0 12 2.1X10 11 3.8xl0 10 1.2xl0 10 5.4xl0 9 3.0xl0 9 

That is, the original estimate [1] 0,(He) * 4 is reasonable, although 
the result does depend logarithmically on xn . 

P 

Although no details were given of the method of calculation, 
Aharonian and Sunyaev [4] have also computed the rate of reaction (1) 
in the high-temperature range. Their lowest T corresponds to the high- 
est value in the above table and the results agree. However the domain 
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0 < 5 is probably most relevant, since it is especially important to 
know when the process begins to be effective as the temperature rises. 
Further details of the calculations described in the present work may be 
found in a separate paper [5]. For example, competing reactions have 
been investigated, such as 

p + 4 He + 20.578 MeV — v 3 He + p + n 

A 3 (11) 

p + 4 He + 19.814 MeV — *- H + p + p 

These are found to contribute only about 1% as much as reaction (1) , due 
essentially to the smaller cross sections for the reactions (11) at low 
energies (just above threshold); this is a phase space effect, 

2 

3. Discussion. In a strong gravitational potential, say, $ ^ 0.1 c , 
it can be expected that a temperature kT 'v* 0.1 Mc^ ^ 100 MeV can be 
attained (M - nucleon mass) . This potential (and temperatures) could 
occur in regions around black holes, for example, and is well beyond 
that necessary to break up He and other nuclei. Matter that has 
attained this temperature can then be expected to be pure hydrogen. How- 
ever, it might not be observed until after it has cooled. We would then 
see a recombination-cascade spectrum with only hydrogen lines after cool- 
ing to^lO^K. While at ^ 10 7 K there would be a pure bremsstrahlung x-ray 
spectrum without x-ray lines from excitations of one- and two-electron 
ions from the K-shell. 
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